Multisite versus multiorbital Coulomb correlations studied within finite-temperature 
exact diagonalization dynamical mean-field theory 
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The influence of short-range Coulomb correlations on the Mott transition in the single-band 
Hubbard model at half-flUing is studied within cellular dynamical mean field theory for square 
and triangular lattices. Finite-temperature exact diagonalization is used to investigate correlations 
within two-, three-, and four-site clusters. Transforming the non-local self-energy from a site basis 
to a molecular orbital basis, we focus on the inter-orbital charge transfer between these cluster 
molecular orbitals in the vicinity of the Mott transition. In all cases studied, the charge transfer 
is found to be small, indicating weak Coulomb induced orbital polarization despite sizable level 
splitting between orbitals. These results demonstrate that all cluster molecular orbitals take part 
in the Mott transition and that the insulating gap opens simultaneously across the entire Fermi 
surface. Thus, at half-filling we do not find orbital-selective Mott transitions, nor a combination 
of band filling and Mott transition in different orbitals. Nevertheless, the approach towards the 
transition differs greatly between cluster orbitals, giving rise to a pronounced momentum variation 
along the Fermi surface, in agreement with previous works. The near absence of Coulomb induced 
orbital polarization in these clusters differs qualitatively from single-site multi-orbital studies of 
several transition metal oxides, where the Mott phase exhibits nearly complete orbital polarization 
as a result of a correlation driven enhancement of the crystal field splitting. The strong single- 
particle coupling among cluster orbitals in the single-band case is identified as the source of this 
difference. 

PACS. 71. 20. Be Transition metals and alloys - 71.27-|-a Strongly correlated electron systems 



I. INTRODUCTION 



Considerable progress has recently been achieved in the 
understanding of the Mott transition in a variety of tran- 
sition metal oxides.^ Whereas density functional theory 
in the local density approximation (LDA) predicts many 
of these materials to be metallic, the explicit treatment 
of local Coulomb interactions via dynamical mean field 
theory (DMFT)^ correctly yields insulating behavior for 
realistic values of the on-site Coulomb energy U. In the 
metallic phase, the non-cubic structure of some of these 
systems gives rise to non-equivalent, partially filled sub- 
bands that are split by a crystal field and exhibit orbital 
dependent electron occupancies. The hallmark of the 
Mott transition of these oxides is that orbital polarization 
can be greatly increased by Coulomb correlations and 
that the insulating phase is nearly completely orbitally 
polarized. For instance, in the case of LaTiOa, the e'g 
bands are pushed above the Fermi level and the remain- 
ing singly occupied ag subband is split into lower and 
upper Hubbard bands^i^ In the case of V2O3, Coulomb 
correlations push the Ug band above the Fermi level, and 
the doubly degenerate e'g subbands exhibit a Mott gap.^'^ 
Also, in the insulating phase of Ca2Ru04, the d^y like 
band is completely filled and the d^z^yz like subbands are 
split into Hubbard bandsi^i^ The common feature of the 
Mott transition in these materials is that the effective 
band degeneracy is reduced from three to two or one, so 
that the critical Coulomb energy is lower than it would be 



if the t2g bands were fully degenerate. On the other hand, 
other materials can exhibit a quite different behavior. 
For instance, orbital polarization in BaVSs was shown 
to decrease with increasing local Coulomb interaction!^ 
Also, the Mott transition in LaVOs and YVO3 occurs 
before orbital polarization is complete.^*' Moreover, in 
a hypothetical tetragonal structure of LaTiOs, relevant 
for heterostructures, the Mott phase is reached when 
n-xz.yz approaches 1/4 and rixy vanisheSf^ Finally, the 
possibility of so-called orbital selective Mott transitions 
in multi-band systems has been discussed extensively in 
the literature ]^'^^i^^'^"^i^^'^^ These different trends under- 
line the remarkably rich physics of Mott transitions in 
multi-orbital materials. 

The aim of this work is to investigate the relationship 
between Coulomb correlations in single-site multi-orbital 
systems as described above to those occurring within 
a single band when inter-site Coulomb correlations 
are taken into account. The influence of short-range 
correlations on the nature of the Mott transition is 
currently of great interest and has been studied by many 

groups ,17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37.38.39.40.41.42 

Here we examine the role of correlation driven orbital 
polarization in the vicinity of the Mott transition. For 
example, it is well known that in a minimal two-site 
cluster model,'*- which permits explicit treatment of 
short-range Coulomb correlations in an isotropic square 
lattice, the Green's function and self-energy become 
diagonal if one transforms the site basis to a diagonal 
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bonding - antibonding molecular orbital basis. In a 
four-site cluster model, diagonality is obtained by trans- 
forming sites to cluster molecular orbitals characterized 
by r = (0,0), X = (7r,0),(0,7r), and M = (7r,7r)4i The 
molecular orbital components of the self-energy provide 
qualitative information on the importance of correlations 
in the corresponding sections of the Brillouin Zone. In 
the case of an isotropic triangular lattice, Green's func- 
tion and self-energy can be diagonalized by an analogous 
transformation to molecular orbitals appropriate for a 
three-site cluster4^ The question then arises whether 
these cluster molecular orbitals in the single-band case 
obey a similar scenario as the multi-orbital systems 
mentioned above. 

Since an approximate momentum variation of the lat- 
tice self-energy in these models can be derived from a 
linear superposition of the respective molecular orbital 
components of the cluster self-energy, the effect of corre- 
lation enhanced orbital polarization is of direct relevance 
for the question of whether the Mott gap opens uniformly 
across the Fermi surface, or whether it opens first in cer- 
tain regions of the Brillouin Zone (e.g. near the so-called 
hot spots) and only at larger U in the remaining regions 
(the so-called cold spots). The latter picture would be 
analogous to the orbital selective Mott transition which 
can occur within single-site DMFT treatments of cer- 
tain multi-band systems.''''"'"'^'"'' iH -J^ii^ Another possibil- 
ity, analogous to multi-orbital materials such as LaTiOa, 
V2O3, and Ca2Ru04, is that a subset of cluster orbitals 
could exhibit a genuine Mott transition, while the re- 
maining ones are pushed above or below the Fermi level 
at about the same critical U . 

To account for inter-site correlations we use DMFT 
combined with finite temperature exact diagonalization 
(ED). 5^ It was recently shown-- that this method can be 
generalized to multi-band materials by computing only 
those excited states of the impurity Hamiltonian that are 
within a narrow range above the ground state, where the 
Boltzmann factor provides the convergence criterion. Ex- 
ploiting the sparseness of the Hamiltonian, these states 
can be computed very efficiently by using the Arnoldi 
algorithm>5i Higher excited states enter via Green's func- 
tions which are evaluated using the Lanczos method. 
This approach has proved to be highly useful for the 
study of strong correlations in several transition metal 
oxides4iiii^i^i^ An important feature of ED/DMFT is 
that low temperatures and large Coulomb energies can 
be reached. The adaptation of single-site multi-orbital 
ED to multi-site single-band systems is discussed in de- 
tail below. In particular, we introduce a mixed site 
- molecular orbital basis which permits a more flexi- 
ble and more accurate projection of the lattice Green's 
function onto the cluster than in a pure site represen- 
tation. Previous multi-site ED/DMFT studies focussed 
on T = o.^i25i37^^^^ rpj^g extension to finite T dis- 
cussed here is especially useful for the evaluation of the 
T/U phase diagram. 

The main result of this work is that in all cluster mod- 



els studied here for half-filled square and triangular lat- 
tices, there is little enhancement of orbital polarization 
in the vicinity of the Mott transition. Thus, despite siz- 
able level splitting between these cluster orbitals, they 
all exhibit Mott gaps at the same critical Coulomb en- 
ergy. As a consequence, the Mott gap in these models 
opens uniformly across the Fermi surface. For the square 
lattice we show explicitly that the Mott gap at the cold 
spot M/2 = (7r/2,7r/2) of the Brillouin Zone is driven 
by Coulomb correlations at the hot spot X = (tt, 0). 
Therefore, there is no orbital selective Mott transition. 
Moreover, there is no evidence for the combination of 
partial band filling and Mott transition in remaining sub- 
bands that is characteristic of single-site DMFT treat- 
ments of the multi-orbital materials LaTiOa, V2O3, and 
Ca2Ru04, as mentioned above. 

The outline of this paper is as follows. Section II dis- 
cusses the theoretical aspects of our cluster ED/DMFT 
implementation of finite temperature exact diagonaliza- 
tion. Section III provides the results for the two-site and 
four-site clusters of the square lattice, and the three-site 
cluster of the isotropic triangular lattice. In Section IV 
we briefly discuss analogies and differences between these 
multi-site correlation effects and those investigated pre- 
viously in single-site DMFT treatments of multi-orbital 
materials. The conclusions are presented in Section V. 



II. MULTI-SITE ED/DMFT 

Let us consider the single-band Hubbard model 

^ = E (^-^J- + ^•^•) + "»T"a (1) 

{ij)a i 

where the sum in the flrst term extends over nearest 
neighbor sites. The hopping integral t will be set equal to 
unity throughout this paper. Thus, the band widths of 
the square and triangular lattices are = 8 and 11^ = 9, 
respectively. Within cellular dynamical mean field theory 
(CDMFT)2^i^ the interacting lattice Green's function in 
the cluster site basis is given by 

(iw„) = ^ + - ^(fc) - S(iLj„)^ (2) 
fc 

where the k sum extends over the reduced Brillouin Zone, 
(jj„ = (2n -I- l)iTkBT are Matsubara frequencies and /j, 
is the chemical potential. The lattice constant is unity. 
t{k) denotes the hopping matrix for the superlattice and 
E(iw„) represents the cluster self-energy matrix. To 
make contact to other recent works (^^iiliiid^dl we con- 
sider here the paramagnetic metal insulator transition. 

In the site basis, the Green's functions for two-site, 
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three-site and four-site clusters have the structure 



^(4) ^ 




(3) 



(4) 



(5) 



with a = Gil, b = G12 and c — G14. Site labels in 
the square lattice refer to 1 = (0,0), 2 = (1,0), 3 = 
(0,1), 4 = (1,1) and in the triangular lattice to 1 = 
(0,0), 2 = (1,0), 3 = (1/2, V3/2). The superscript de- 
notes the cluster size ric in the square lattice {ric — 2 
or Tic = 4) or triangular lattice {ric = 3), respectively. 
In the site bases, the corresponding self-energy matri- 
ces S^""^^ (i^^n) have the same symmetry properties as the 
Green's functions. 

A key aspect of DMFT is that, to avoid double- 
counting of Coulomb interactions in the quantum impu- 
rity calculation, it is necessary to remove the self-energy 
from the cluster in which correlations are treated exactly. 
This removal yields the Green's function 



Go(iw„) = [G{iu!n) ^ + S(iw„)] 



(6) 



These matrices also exhibit the symmetry properties 
specified above. 

For the purpose of the ED calculations it is convenient 
to transform the site bases into molecular orbital bases 
in which the Green's functions and self-energies become 
diagonal. For the two-site cluster, molecular orbitals are 
given by the bonding - anti-bonding combinations 01^2 — 
(|1) ± |2))/V2. For the four-site cluster, they are formed 
by the plaquettes 0i = (|1) + |2) + |3) + |4))/2, 02 = 
(|1) + |2) - |3) - |4))/2, 03 = (|1) - |2) + |3) - |4))/2, 
04 = (|1) - |2) - |3) + |4))/2. Finally, for the three-site 
cluster of the triangular lattice they can be written as: 
01 = (|1) + |2) + |3))/V3, 02 = (-2|1) + |2) + |3))/^/6, 
03 = (1 2) — \3))/\/2. In theses cluster molecular orbital 
bases, the above Green's functions take the form 



0.3 



to 0.2 



0.1 



0.3 



0.2 



0.1 



0.3 



0.2 



0.1 



1 1 1 1 1 
square lattice 11 

2 sites / \ 


f Pi 

J\---'' 1 1 1 


iV 


-4-2 2 

(0 


4 


1 1 j 1 
square lattice 11 




4 sites / \ 




^■"■^ / \ P2,3 




f P1 ^'',>--- ><" 


^ P4 1 



triangular lattice 




3 sites 






•'P2,3'l 


P- 



-2 

(0 



FIG. 1: (Color online) Total density of states p(cj) and molec- 
ular orbital components pm(i^) for two-site and four-site clus- 
ters of square lattice (top panels) , and of three-site cluster of 
triangular lattice (bottom panel). For clarity, the molecular 
orbital components are divided by ric. 



In the site basis, the local density of states in the non- 
interacting limit is given by 



G(2) = 
G(3) = 

G(4) ^ 
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(7) 

(8) 

(9) 



The self-energies S(ia;„) and Green's functions Go(ia;„) 
can be diagonalized in the same fashion. We denote these 
elements as Gm(i^n), Sm(iw„) and Go^mii^n)- 



Pii{uj) = --lmGii(w) 

TT 



(10) 



with E = 0. Since we consider isotropic clusters, all 
sites are equivalent, so that pii{uj) coincides with the 
density of states p{uj). In the molecular orbital basis, 
the density of states components Pm{oj) have different 
shapes and different centroids, analogous to the crystal 
field split density of states components of many transition 
metal oxides. Figure 1 shows these densities for the two- 
site and four-site clusters of the square lattice and the 
three-site cluster of the triangular lattice, as described 
above. According to Eqs. ([7][n]), the cluster molecular 
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orbital densities of states are given by 
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where p|j''' are the site components for cluster Uc. 

From these cluster molecular orbital densities of states 
an approximate momentum variation across the Brillouin 
Zone can be constructed (see below) . For instance, in the 
case of the square lattice with ric — 4, densities associ- 
ated with the high-symmetry points of the original lattice 
are given by prito) = piiuj), Px{t^) = P2(t^) = P3(w) and 
Pm{oj) = Piil^)- At M/2 = (7r/2,7r/2), the density of 
states corresponds to the local density p{lo) = pii[tjj) = 
{pr + Pm + 2px)/4:. Note, however, that all molecu- 
lar orbital densities extend across the entire band width. 
Thus, they are not identical with those sections of the lo- 
cal density of states that originate in momentum regions 
surrounding the high-symmetry points, as would be the 
case in the dynamical cluster approximation (DCA)^'^ 

We now project the Green's function Go{iuJn) defined 
in Eq. ^ onto a cluster consisting of Uc impurity levels 
and Ub bath levels. The total number of levels is Us — 
nc + ni,. In the site basis we have 



Go{iLu„) ~ Go(iw„) 



(14) 



where h is the non-interacting impurity cluster Hamilto- 
nian and T{iLUn) the hybridization matrix describing the 
coupling between impurity cluster and bath. Thus, 




(15) 
(16) 

(17) 



For the square lattice we choose eo = and for the trian- 
gular lattice eo — —0.83, so that the Fermi level coincides 
with u> — 0. 

Instead of expressing the non-diagonal hybridization 
matrix T{iujn) in a site basis, it is convenient to go over to 
the molecular orbital basis in which G'o(iwn) is diagonal. 
Assuming that each component G'o.m(ia;„) couples only 



with its own bath, we have 
Go,m(iw„) « Go'„(ia;„) 

= \iuj.n + p- 



(18) 



where em represents an impurity level, efe the bath levels, 
and Vmk the hybridization matrix elements. The incor- 
poration of the impurity level e,„ ensures a much better 
fit of Go^m{i^n) than by projecting only onto bath or- 
bitals. 

For instance, for ric — 3 and ~ 12 (i.e., three bath 
levels per impurity orbital), each component Go,m(iw„) 
is fitted using seven parameters: one impurity level £„, 
three bath levels and three hopping integrals Vmk- 
Since according to Eq. ([5]) there are two independent 
functions, we use a total of 14 fit parameters to repre- 
sent these two Go,m components. This procedure allows 
for a considerably more flexible projection of the Green's 
function matrix Go(jWn) onto the bath. In a site ba- 
sis for an isotropic triangular lattice (taking again three 
bath levels per site) one would have instead only 6 fit 
parameters if each site couples to its own bath. Effec- 
tively, therefore, the molecular orbital basis accounts for 
several additional cross hybridization terms as well as 
internal cluster couplings (see below). Moreover, it is 
much more reliable to fit the two independent molecular 
orbital components Go,m(i^n) than a non-diagonal site 
matrix Go,ij(iti;„) with an equivalent number of parame- 
ters. Analogous considerations hold for the two-site and 
four-site clusters of the square lattice. For example, for 
Tie — 4:, Hs — 12 there are two independent functions Go,i 
and Go, 2 (Go, 4 is related to Go,i), giving a total of 10 fit 
parameters, compared to only two parameters in a simple 
site picture with fourforld and particle hole symmetry. 

Figure 2 illustrates the typical quality of the projec- 
tion of the lattice components of Go,m(«^^n) onto the 
bath for ric = 3 and rig = 9. Thus, although only 
two bath levels per orbital are included (i.e., using five 
parameters per orbital), the fit of both real and imag- 
inary parts is excellent. For these cluster sizes and 
low temperatures, iterations take only a few minutes. 
Fits of similar quality are achieved for multi-orbital 
materials j^iiiii^i^'^ To achieve even better agreement 
at low frequencies, it is preferable to minimize not the 
bare difference Go^mi'i^n) — Gq m{iujn) but to divide these 
functions first by ujn- 

We now discuss the evaluation of the finite tempera- 
ture interacting cluster Green's function. If this step is 
carried out in the diagonal molecular orbital basis, the 
Coulomb interaction must be expressed as a matrix con- 
taining many inter-orbital components. For ric = 4 it 
can be easily shown that Umim2m3mi = C//4 for 64 of the 
possible 256 configurations. All other matrix elements 
vanish. This step can be circumvented by working in a 
mixed basis consisting of cluster sites i and bath orbitals 
k. We illustrate this procedure here for the triangular 
lattice with nc = 3. Let us denote the transformation 
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FIG. 2: (Color online) Projection of lattice Green's function 
components Go,m{ii^n) onto bath for — 3, U — 9, T = 
0.02. Upper panel: ImGo,™, lower panel: ReGo,m. Red 
curves: lattice Green's functions, blue curves: approximate 
expression, right-hand side of Eq. (fTSjl . 



between sites and orbitals as T^"''\ where 



T, 



(3) _ 



1/^3 -2/V6 
1/V3 1/a/6 I/V2 
1/V3 l/\/6 -l/\/2. 



(19) 



In this mixed basis, the effective site block of the cluster 
Hamiltonian becomes 



e' t' t' 
t' e' t' 
t' t' e' 



(20) 



with e' = eo + (ei + 2e2)/3 and t' = t + {ti - e2)/3, 
where eo and t are the elements of the original cluster 
Hamiltonian defined in Eq. (|16p . The new terms involv- 
ing the molecular orbital cluster levels arise from the 
projection specified in Eq. p8|) . In the mixed basis, the 
hybridization matrix elements Vmk between cluster and 
bath orbitals introduced in Eq. are transformed to 
new hybridization matrix elements between cluster sites 
i and bath orbitals k. They are given by 



mk • 



(21) 



Using the elements e', t' and V-/, together with the 
on-site Coulomb energy U, the non-diagonal interacting 



cluster Green's function at finite temperature is derived 
from the expressionSiS 



G^-(zc.„) 



z ^ 



E,, - E„ + iw„ 



(22) 



where E^, and denote the eigenvalues and eigen- 
vectors of the impurity Hamiltonian, f3 — l/ksT and 
Z = exp(— /3i5^) is the partition fvmction. At low 
temperatures only a relatively small number of excited 
states in few spin sectors contributes to Gfy They can 
be efficiently evaluated using the Arnoldi algorithm.— 
The excited state Green's functions are computed using 
the Lanczos procedure. Further details can be found in 
Refi^. The non-diagonal elements of Gfj are derived by 
first evaluating the diagonal components C?^,- and then 
using the relation^ 



Since G^j ~ G^-, this yields: 



Gf + G 



31 



G 



33 ■ 



GC 



33 



(23) 



(24) 



For the two-site cluster, we have used 3 or 4 bath levels 
for each impurity orbital [ug = 8 or 10), for the three-site 
cluster 2 or 3 bath levels per impurity orbital (n^ = 9 or 
12), and for the four-site cluster 2 bath levels per impu- 
rity orbital {rig — 12). Gfj{iujn) obeys the same sym- 
metry properties as the lattice Green's functions given 
in Eqs. (l3]-[5]). It therefore can be diagonalizcd as indi- 
cated in Eqs. ([7H51). We denote these diagonal elements 
as G^(zci;„). For Uc = 4, we have checked that the evalu- 
ation of G^(iw„) in the non-diagonal site - orbital basis 
and in the diagonal molecular orbital basis yield identical 
results. 

The key assumption in DMFT is now that the result- 
ing impurity cluster self-energy is a physically reasonable 
representation of the lattice self-energy. Thus, using a re- 
lation analogous to Eq. ([5]), we find 



l/G^oU'uj„ 

Sm(iw„). 



(25) 



After transforming I]m(iw„) back to the non-diagonal site 
basis, it is used as input in the lattice Green's function 
Eq. ^ in the next iteration step. 

To summarize the procedure discussed above, the 
multi-site ED/DMFT calculation consists of the follow- 
ing steps: 

(a) evaluate the lattice Green's function Gij{iLUn), 
Eq. ([2]), in the non-diagonal site basis, using as input 
the self-energy obtained in a previous iteration step. 
The entire iteration procedure is started at small U with 
S = 0. 
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(b) transform Gy and to the diagonal molecular 
orbital basis and compute the components Go,m- 

(c) project the Go,m(ici^n) onto independent baths to 
determine e„i, Cfe and Vmk as indicated in Eq. (HH]). 

(d) from the fit parameters and Vmk determine the 
Hamiltonian matrix elements e', t' and V^f. in the mixed 
site - orbital basis. 

(e) evaluate the non-diagonal cluster Green's function 
Gfj{iujn) using the Arnoldi and Lanczos methods. 

(f) transform this Green's function to the diagonal 
orbital basis and compute the cluster self-energy compo- 
nents E^(icj„) defined in Eq. (^5)1 . 

We emphasize that ED/DMFT involves, at each iter- 
ation, two projections: (1) The lattice Green's function 
Go is projected onto the cluster Green's function Gq, as 
indicated in Eq. ^TE\i . By definition, Go has a continuous 
spectrum at real frequencies, while Gq is discrete. (2) 
The cluster self-energy E'^', which evidently has a dis- 
crete spectrum at real to, is used as an approximation of 
the lattice self-energy S, which by definition is contin- 
uous along the real frequency axis. Thus, both projec- 
tions. Go « Gq and E'^' « E, rely on the well-known fact 
that continuous and discrete spectra at real to can yield 
nearly identical distributions at Matsubara frequencies. 
Since the cluster size determines the number of discrete 
spectral features of Gq and E^' , there exists evidently an 
infinite number of discrete spectra which may in principle 
be used to represent the continuous spectra of the lattice 
quantities Go and E. 



III. RESULTS AND DISCUSSION 
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Figure 3 shows the occupancies of the cluster molecular 
orbitals for three cluster sizes as functions of increasing 
U. For the square lattice the Mott transition occurs near 
Uc2 ~ 5.5, while for the triangular lattice Uc2 ~ 9.5 (see 
below). Evidently all orbital occupancies vary smoothly 
across the transition. There is no indication of orbital se- 
lective Mott transitions, nor for complete filling or emp- 
tying of any orbitals at large Coulomb energies. This 
behavior differs qualitatively from the one found in ma- 
terials such as LaTiOa, V2O3, and Ca2Ru04, where the 
insulating phase exhibits nearly complete orbital polar- 
ization as a result of a Coulomb driven enhancement 
of the crystal field splitting between t2g orbit als^i^i^i^!^ 
(see following section). 

The fact that all cluster orbitals remain partially oc- 
cupied across the transition implies that the gap opens 
simultaneously in all orbitals. This can be seen most 
clearly in the spectral distributions, as shown in Figure 4 
for the square lattice with ric = 4. Since we are here con- 
cerned with the transition from metallic to insulating be- 
havior, we show the spectra obtained from the interacting 
cluster Green's function, Ay(w) = — (I/tt) ImG^j(w-|-i(5), 
Eq. (im, with S = 0.1. These spectra can be evalu- 
ated without requiring analytic continuation from Mat- 



FIG. 3: (Color online) Occupancies of cluster molecular or- 
bitals (per spin) for two- and four-site clusters of square lat- 
tice (top panels) , and of three-site cluster of triangular lattice 
(bottom panel) as functions of Coulomb energy U. The tem- 
perature is r = 0.02. The sum of these molecular occupancies 
is nc/2 and their average is 0.5. The Mott transition for in- 
creasing U occurs at Uc2 ~ 5.5 for the square lattice (ric = 2 
and Tic — 4), and near Uc2 ~ 9.5 for the triangular lattice 
(jic = 3). 



subara to real frequencies. Using the transformations 
indicated in Eqs. ([9]) and ((13]), the cluster molecular or- 
bital densities are given by Ai — An + 2Ai2 + A14, 
A2 = A3 = All - Ai4 and A4 = An - 2Ai2 + A14. The 
total density ^(w) = Am(w)/4 coincides with the on- 
site distribution Aii{lo), which also represents the den- 
sity corresponding to M/2 = (7r/2,7r/2) (see below). To 
our knowledge, this orbital decomposition has not been 
addressed before. Clearly, all orbitals contribute to the 
spectral weight at Ep in the metallic phase, as well as 
to the lower and upper Hubbard bands in the insulat- 
ing phase. The spectral distributions shown in the upper 
panel are consistent with those by Kyung et ali^ and 
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FIG. 4: (Color online) Upper panel: total spectral distri- 
bution below and above the Mott transition (f/c « 5.5) for 
the square lattice with nc = 4 at temperature T = 0.02: 
U = 5 (red curve) and U — 6 (blue curve). Lower two panels: 
molecular orbital contributions at (7 = 5 and U — G; blue 
curves: ^1(0;) (associated with F), green curves: A4{uj) (M), 
red curves: A2,3{io) (X) (broadening 5 = 0.1). 



Zhang and Imada^ within ED/DMFT at T = 0. The 
spectra reveal a characteristic four-peak structure, con- 
sisting of low- frequency peaks limiting the pseudogap due 
to short-range correlations and high-frequency peaks as- 
sociated with the Hubbard bandsi^^ii^i^ The small peak 
at Ep in the metallic phase at J7 = 5 appears only at 
finite T. It vanishes at T = 0. At such low frequencies, 
however, ED finite-size effects cannot be ruled out. 

Similar results are obtained for the square lattice in the 
two-site cluster model. Results for the triangular lattice 
with Tic — 3 are shown in Figure 5. According to Eqs. ([8]) 
and (|12p the cluster molecular orbital densities are given 
by Ai = All + 2^12 and A2 ^ A3 ^ An - A12. As for 
the square lattice, the Mott gap opens simultaneously 



FIG. 5: (Color online) Spectral distributions of cluster molec- 
ular orbitals below and above the Mott transition {Uc ~ 9.5) 
for the triangular lattice with nc = 3 at temperature T = 0.02. 
Upper panel: U = 9.0; lower panel: U — 9.6. Blue curves: 
j4i(w), red curves: ^2,3(0;). 



in all orbitals at about the same critical U, and all or- 
bitals contribute to the lower and upper Hubbard bands. 
While the metallic phase of the unfrustrated square lat- 
tice close to the transition exhibits a pseudogap due to 
short-range antiferromagnetic correlations r^dldidS this 
phenomenon is absent in the triangular lattice as a result 
of geometrical frustrationi^ 

To illustrate the first-order nature of the metal insula- 
tor transition we show in the upper panel of Figure 6 the 
spectral weight of the ric = 3 cluster orbitals at Ep — 
as a function of U. The lower panel shows the average 
double occupancy docc = J2mi''^m'\nmi) /3- Both quan- 
tities exhibit hysteresis for increasing and decreasing U, 
indicating the coexistence of metallic and insulating so- 
lutions. The complete T/U phase diagram will be pub- 
lished elsewhere. The phase diagram for the isotropic 
square lattice with nc = 4 was recently mapped out in 
detail by Park et al.''^ 

The cluster molecular orbital components of the self- 
energy and Green's function may be used to derive an 
approximate expression for the momentum variation of 
the lattice self-energy and Green's function in the original 
Brillouin Zone>22i 

I](fc,zc^„) « 1 Ve^^-(^'-^^)S,,(zc^„) (26) 



where Ri are the cluster site positions and T^ij are the site 
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the diagonal on-site element of the cluster self-energy, 
which is identical with the local lattice self-energy. Thus: 
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FIG. 6: (Color online) Hysteresis behavior of spectral weights 
A,n{0) of cluster molecular orbitals at Ef ~ (upper panel) 
and average double occupancy (lower panel) as functions of 
Coulomb energy for triangular lattice with tIc = 3 (T = 0.02). 
Red curves: increasing U, blue curves: decreasing U. 



components of the self-energy. An analogous expression 
holds for the lattice Green's function. For the clusters 
discussed above this superposition implies that the Mott 
gap opens uniformly along the Fermi surface since all 
orbitals undergo a common transition. 

Writing the site elements in terms of the orbital com- 
ponents, one has for ric — 4: E(r, icij„) = Ei(icij„), 
E(X, ia;„) = S2(iw„) = S3(iw„) and S(M, jw„) = 
1^4(iLUn)- In agreement with results of previous 
authorsHiSS^di^^^ behavior of S(r, icj„) 

and S(M, iw„) near the Mott transition to differ qualita- 
tively from that of iujn) (not shown here): Whereas 
ImE(X, iwn) exhibits ~ w„ variation at low frequencies 
in the metallic phase and 1 /uj„ variation in the insulat- 
ing phase (the real part vanishes because of particle hole 
symmetry, see Fig. 1), Im S(r, and Im S(M, iujn) re- 
main ~ LJn in both phases, but their real parts increase 
rapidly across the transition. 

Although this behavior might suggest a Mott transi- 
tion for the X cluster orbitals combined with a band- 
filling or band-emptying mechanism for the F and M or- 
bitals, the orbital occupancies (Fig. 3) and spectral dis- 
tributions (Fig. 4) demonstrate this not to be the case. 
Moreover, for the crucial question of whether or not the 
Mott gap opens simultaneously across the Fermi surface, 
it is important to compare the self-energy at X with its 
behavior at M/2, where the electron band also crosses 
Ep- According to Eq. (pS)) . S(M/2,ia;„) coincides with 



T,{M/2,iuj„ 



1 



[I](r, iuJn) + S(Af, iuJn) 
+ 2S(X,iw„)]. 



(27) 



The real parts of the first terms on the rhs cancel since 
the corresponding density of states components are mir- 
rors of each other (see Fig. 1). Thus, S(M/2,zu;„) 
and T,{X,iujn) are purely imaginary because of particle 
hole symmetry. The above relation demonstrates that 
Coulomb correlations at the cold spot A//2 are essentially 
driven by those at the hot spot X. In fact, the magnitude 
of the self-energy at M/2 is a factor of 2 smaller than the 
singular term at X, with weak additional, non-singular 
contributions associated with F and M. 

Along the Fermi surface between X and M/2 (i.e. for 
kx + ky = tt) the self-energy is given by 

T,{k,iujn) — cos^(fcj;)I](A", iw„) + sin^(fcj;)S(A//2, ia;„). 

(28) 

This function is imaginary, i.e., there are no band shifts 
due to a finite real part of the self-energy. Thus, within 
the four-site cluster DMFT, the opening of the Mott gap 
on the entire Fermi surface is determined solely by the 
singularity of Im I](X, i(jj„). 

The above analysis leads to a surprisingly simple pic- 
ture for the momentum variation of correlations along 
the Fermi surface. It consists of two sinusoidal contri- 
butions: The singular X term oscillates with amplitude 
1 at X and 1/2 at M/2, and the non-singular term due 
to the F, M orbitals oscillates with amplitude 1 at M/2 
and zero at X. 

This is illustrated in Figure 7 which shows the lattice 
self-energy and Green's function at X and M/2 below 
and above the Mott transition. In the metallic phase, 
the self-energy at X is seen to be larger than the one 
at M/2. Thus, quasi-particle lifetimes decrease between 
M/2 and X. Nevertheless, the singular behavior at M/2 
is governed by the one at X. Accordingly, the Green's 
function at X and AI/2 displays a change from metallic 
to insulating behavior at the same Coulomb energy. 

To understand the strong coupling between different 
sections of the Brillouin Zone it is important to recall 
that, in the diagonal cluster molecular orbital basis, the 
single-particle part of the lattice Hamiltonian appearing 
in Eq. ^ is not diagonal. Thus, the orbital component 
of the lattice Green's function Gra{i^n) is infiucnccd not 
only by the corresponding self-energy Em(zw„), but by 
the other orbital elements as well. This point becomes 
clear if we compare Gm{i^n) with the approximation 



Gni{i(^n) ~ duj 



+ /i - W - Y^raiioJu) ' 



(29) 



Figure 8 shows that in the metallic region aXU — b there 
is little difference with regard to the actual Gmii^n)- 
Also, at [/ = 6 the key components responsible for the 
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FIG. 7: (Color online) Lattice self-energy, Eq. (I26|l . (upper 
panel) and Green's function (lower panel) ai X — (tt, 0) (red 
curves) and M/2 — (7r/2, 7r/2) (blue curves) for metallic phase 
{U = 5) and insulation phase {U = 6) of square lattice {ric — 
4) at T = 0.02. Both iuj„) and G{k,iuJn) are purely 
imaginary between X and M/2. The sinusoidal variation of 
the self-energy between these points is given by Eq. (|28|) . 



metal insulator transition, namely G2,z{ii^n) correspond- 
ing to X = (tTjO), are well represented by this approxi- 
mation. The T,M components, Gi,4(ict;„), however, do 
not reveal insulating behavior since at small cj„ the imag- 
inary parts do not extrapolate to zero. This demonstrates 
that the Mott gaps seen in Ai. 4(0-1) in Figs. 4 and 5 are 
not caused by the rapidly varying real parts of Y,i^i{iuJn). 
Instead, the gaps at T and M are driven by the singular 
behavior of Im S2,3(io;„) which contributes to Ai^^^iio) via 

the non-diagonal elements of t{k).^^ 

If the approximate components obtained via Eq. (j29p 
were used to generate the spectral distributions Am{uj), 
it is clear that only A2,3(w) corresponding to X would 
exhibit the Mott transition, while ^1^4(0;) would retain 
considerable metallicity. This implies that the physics at 
the cold spot M/2 is incorrectly represented via Eq. (p9)) . 
suggesting that the gap at M/2 does not open at the same 
U as at X. Instead, as argued above, the metal insulator 
transition at M/2 is caused by the same self-energy terms 
as at X, i.e., the Mott gap opens uniformly. 

As pointed out earlier, in the molecular orbital ba- 
sis the Coulomb matrix has a large number of non-zero 
elements. Thus, in this basis there is not only single- 
particle hybridization arising from t(fc), but also strong 
inter-orbital Coulomb repulsion. Nevertheless, the above 



FIG. 8; (Color online) Lattice Green's function Eq. (2) in 
He = A cluster molecular orbital basis for f/ = 5 (upper panel) 
and ?7 = 6 (lower panel) at T = 0.02. Red curves: diago- 
nal version of Eq. ([21; blue curves: approximate expression, 
Eq. ((291). 



analysis reveals that, although these Coulomb interac- 
tion terms are properly taken into account, the spectral 
distributions at P, M and M/2 do not exhibit a Mott gap 
unless the non-diagonal terms of t{k) in the orbital ba- 
sis are included. It is therefore the single-particle part of 
the Hamiltonian that provides the correct connection be- 
tween the self-energy components and thereby generates 
the true momentum variation of the spectral distribution 
of the single band. 

Because of the strong coupling between orbitals, the 
notion that some of these orbitals undergo a Mott tran- 
sition while others do not, does not appear appropriate. 
As shown consistently by all orbital-resolved spectra in 
the present multi-site ED /DMFT study, there is a single 
Mott transition common to all cluster molecular orbitals, 
implying a simultaneous opening of the Mott gap across 
the entire Fermi-surface. Nevertheless, in agreement with 
previous authors, we find pronounced momentum varia- 
tion of quasi-particle properties close to the Mott transi- 
tion. 

In view of the approximate nature of the momentum 
variation of the lattice self-energy and Green's function 
derived within the CDMFT it would be very interest- 
ing to compare the above results with analogous ones 
obtained within the DCA)2£i2& in particular, since the 
cluster molecular orbital components of the density of 
states, as stated above, differ appreciably between these 
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FIG. 9: (Color online) Subband occupancies obtained within 
finite temperature multi-band ED/DMFT as functions of U. 
Upper panels: LaTiOa {3d^ configuration) for full Hund ex- 
change with J = 0.65 eV— ; V2O3 configuration) for full 
Hund exchange (blue dots) and Ising-like exchange (red solid 
dots) with J — 0.7 eV— ; red empty dots: quantum Monte 
Carlo DMFT results.^ Lower panel: Ca2Ru04 (4d"' configu- 
ration) for full Hund exchange with J = 17/4 and A = 0.4 eV 
t2g crystal field splitting.— 



two cluster DMFT schemes.~ This comparison will be 
addressed in a future publication. 



The scenario discussed in this section differs strikingly 
from the one found for several multi-orbital materials 
which have been studied previously by various groups. 
For the sake of comparison we review some of these sys- 
tems in the following section. 



FIG. 10: (Color online) Schematic illustration of correlation 
driven enhancement of orbital polarization. Upper row: crys- 
tal field split t2g LDA densities of states for occupancies n = 1, 
2, and 4, corresponding to LaTiOa, V2O3 and Ca2Ru04, 
respectively. Blue curves: singly-degenerate ag band, red 
curves: doubly-degenerate e'g bands. In the case of Ca2Ru04, 
a refers to dxy, e to dxz,yz- The vertical bars denote the Fermi 
level. Lower row: orbitally polarized Mott phase, n = 1: 
empty e'g bands, lower and upper Hubbard peaks of half-filled 
ag band; n = 2: empty ag band, lower and upper Hubbard 
peaks of half-filled e'g bands; n = 4: filled ag band, lower and 
upper Hubbard peaks of half-filled e'g bands. 



IV. COMPARISON WITH MULTI-ORBITAL 
SYSTEMS 

During the recent years single-site multi-orbital DMFT 
has been used extensively to investigate the metal insula- 
tor transition of a variety of materialsii Here, we briefly 
discuss some of these systems which are regarded as typ- 
ical Mott insulators, and which all exhibit characteristic 
changes of the electronic structure as the metallic phase 
is replaced by the insulator at large Coulomb energies. 

Figure 9 shows the correlation driven enhancement of 
orbital polarization for LaTiOa, V2O3 and Ca2Ru04. 
Because of the orthorhombic structure of LaTiOa, LDA 
calculations reveal that the ag subbands of the t2g sector 
are slightly more occupied than the two e'g componentsj^ 
Local Coulomb interactions enhance this t2g crystal field 
splitting, so that at the Mott transition close ioU = 5 cV 
the e'g bands become nearly empty and the ag band half- 
filled.''-"' In the case of V2O3, the corundum lattice struc- 
ture ensures that the doubly degenerate e'g bands have 
slightly larger binding energy than the ag bands. ^ With 
increasing Coulomb interaction this crystal field splitting 
is strongly enhanced, until in the range [/ w 5 . . . 6 eV, 
the ttg bands are pushed above the Fermi level and the 
e'g bands become half-filled.^1^1^ Finally, in the case of 
Ca a;Sr2_£!;Ru04, Sr substitution via the smaller Ca ions 
gives rise to an enlarged crystal field splitting between dxy 
and dxz,yz like subbands.^'' Coulomb correlations increase 
this splitting, until in the Mott phase the dxy bands are 
fully occupied and the half-filled dxz,yz bands are split 
into lower and upper Hubbard bands4 

Schematically, the uncorrelated densities of states 
of these transition metal oxides and the spectra de- 
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rived within single-site multi-orbital DMFT for realistic 
Coulomb energies are shown in Figure 10. Note that in 
the metallic phase, both orbital symmetries contribute 
to the spectral weight at the Fermi level. In the Mott 
phase, the gap involes transitions between states of op- 
posite symmetry character. 

Despite the different subband occupancies of these ma- 
terials, they exhibit a similar correlation driven enhance- 
ment of orbital polarizationi^ By pushing some sub- 
bands above or below the Fermi level, the effective degen- 
eracy is reduced from three to two or one. The Mott tran- 
sition therefore occurs at lower critical Coulomb energy 
than in a cubic environment with equivalent subbands. 

According to Figure 1, the cluster molecular orbital 
densities of the single band Hubbard model in a multi- 
site picture also exhibit a substantial splitting relative 
to the total band width. Nevertheless, Coulomb correla- 
tions in these cases only lead to moderate charge trans- 
fer between these cluster orbitals, as demonstrated by 
the results given in Fig. 2. The main physical reason 
for this qualitative difference with respect to the multi- 
orbital materials is the strong single-particle hybridiza- 
tion among orbitals so that any tendency towards orbital 
selective Mott transitions is suppressed. For the same 
reason, partial band filling or emptying, with a Mott 
transition in the remaining subset of bands, as found in 
several multi-orbital materials, is also absent. 

A certain amount of intcr-orbital hybridization ex- 
ists also in multi-orbital systems since the single-electron 
Hamiltonian in the orbital basis is not diagonal through- 
out the Brillouin Zone. This residual coupling, however, 
is much weaker than in the single-band multi-site system, 
so that Coulomb correlations can indeed lead to nearly 
complete orbital polarization. This is supported by the 
DMFT results for LaTiOs and V2O3 which have been 
studied both by evaluating the lattice Green's function 
via Eq. ^ (see Refs^i^) and via the approximate ver- 
sion, Eq. (1^51) (see Refsi^i^i^). In these systems both 
formulations give very similar results, in particular, both 
confirm the scenario of strong orbital polarization. 



V. CONCLUSION 

Cellular DMFT combined with finite temperature ex- 
act diagonalization has been used to investigate the in- 



fluence of short range correlations on the Mott transition 
in the single-band Hubbard model. Both square and tri- 
angular lattices at half-filling were studied. A mixed ba- 
sis consisting of cluster sites and bath molecular orbitals 
was shown to provide an efficient method for the evalu- 
ation of the cluster self-energies and Green's functions. 
Since in the cluster molecular orbital representation these 
quantities become diagonal, an intriguing analogy exists 
between Coulomb correlations in these multi-site single- 
band systems and several multi-orbital materials which 
were studied previously within single-site DMFT. 

In remarkable contrast to LaTiOa, V2O3, and 
Ca2Ru04, which exhibit pronounced orbital polarization 
at the Mott transition, the single-band systems show very 
little correlation driven enhancement of orbital polariza- 
tion. Thus, all cluster molecular orbitals take part in the 
metal insulator transition. Moreover, the transition oc- 
curs at the same critical U for all cluster orbitals. Since 
an approximate momentum variation of the lattice self- 
energy and Green's function can be constructed from 
a superposition of these molecular orbital components, 
this finding yields the important result that the Mott 
gap opens simultaneously across the entire Fermi surface. 
Thus, for both square and triangular lattices at half fill- 
ing, there is no orbital selective Mott transition, where 
certain sections of the Brillouin Zone would open a gap at 
lower Coulomb energy than other parts. Moreover, there 
is no evidence for the combination of subband-fiUing and 
Mott transition in other subbands, that is characteristic 
of the multi-orbital materials mentioned above. It would 
be of great interest to investigate whether these findings 
also hold at a finer momentum resolution which would 
require cluster sizes larger than Uc = 4. 
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